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Abstract 

We apply the technique of quasi-adiabatic continuation to study systems with continuous sym- 
metries. We first derive a general form of Goldstone's theorem applicable to gapped nonrelativistic 
systems with continuous symmetries. We then show that for a fermionic system with a spin gap, 
it is possible to insert 7r-flux into a cylinder with only exponentially small change in the energy of 
the system, a scenario which covers several physically interesting cases such as an s-wave super- 
conductor or a resonating valence bond state. 
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Systems with a spin gap are believed to have a number of interesting properties. One 
expects these systems to have short-range spin correlations, but in return exotic properties 
such as fractionalization and topological order often appear In this paper, we apply 
the recently developed technique of quasi- adiabatic continuation [2] to study these systems, 
generalizing ideas developed in [3j]. We start by proving a version of Goldstone's theorem, 
demonstrating exponential decay of spin correlations in a gap ped system with continuous 
symmetries. This result is stronger than previous results % |5|, |6( which proved exponential 
decay of connected correlation functions in gapped systems, as in this case we also show that 
for a system with multiple ground states the average of the spin correlations over different 
ground states is exponentially small. 

We then consider the problem of flux quantization in electron systems with a spin gap 
(without requiring that there be any gap to spinless excitations). We consider a system 
of spinful fermions which is periodic in at least one direction, and show that if the system 
has a spin gap then there is a state which has an expectation value of the energy for a 
Hamiltonian with ir flux which is close to the ground state energy of the Hamiltonian with 
no fluxj?}]. This result covers various interesting physical cases such as superconductors s| 
and resonating valence bond states^]. 

Throughout, we consider local Hamiltonians on a lattice as defined below. We begin 
by defining some notation and defining conditions for a Hamiltonian to be local. Then we 
review the technique of quasi-adiabatic continuation, and give a brief discussion of Berry's 
phase for this continuation. We then consider Goldstone's theorem and the flux periodicity 
problem. 

We use labels i,j, ... to denote lattice sites, and introduce a metric dist(z, j). We assume 
that there is a finite dimensional Hilbert space on each site. An example of a local Hamil- 
tonian is a Hamiltonian of the form 7~t = J2iHi, where each Hi has bounded operator norm: 
< J for some J, and where each Hi has support on a set of sites within some distance 
R of site i. This includes lattice spin and fermion models with finite range interactions. 

In general, we define a Hamiltonian to be local as follows: we write the Hamiltonian as 
TC = J2z hz where Z are different sets of sites on the lattice and hz is a hermitian operator 
with support on set Z. For a Hamiltonian to be local, we require that there exist constants 
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(jl, Si such that for all i, 



\\hz\\ \Z\ exp[/zdiam(Z)] < s 1 < oo, (1) 



where \Z\ is the total number of sites within set Z (the cardinality). For Goldstone's theorem 
we also impose an additional finite dimensionality requirement, Eq. (T2TT) . discussed below. 



Eq. (pQ) implies a Lieb- Robinson bound 5|, |6|, llCj : there exists a velocity t> and constant c 



such that for any two operators A, B defined with support on sets X, Y, 

\\[A(t),B]\\ <cx |X|||A||||S||exp[-/idist(X,r)](exp[-y//|t|] - 1), (2) 

where A(t) = exp[i7it]A exp[— iHt]. This bound is useful for \t\ < l/v, and implies that for 
any // < /i, we have 

|p(t),£]|| <cx |X| pim^ll exp[- y u / (dist(X,F)], (3) 
so long as t < l/v', with 

i/ = vn/(n-n'). (4) 

We will use the bound in the form ([3]) throughout. 

Finally, we emphasize that our results are derived for finite volume systems. The bounds 
that we find, however, for Goldstone's theorem, do not depend on volume, while we will 
explicitly discuss the volume dependence for the flux quantization result. Stronger Goldstone 
bounds may be possible in some cases in the infinite volume limit, as there exist some 
systems for which, for example, there is a degenerate ground state for a finite size system 



but a unique ground state in the infinite system|ll[]; we do not consider these possibilities 
further. There also exist cases in which there are two approximately degenerate low energy 
states, which become exactly degenerate in the thermodynamic limit, with a gap to the rest 
of the spectrum 1 121]; our results are applicable to this case, as in the proof of the Goldstone 
theorem we do not require exact degeneracy of the low-lying states. 

I. QUASI- ADIABATIC CONTINUATION 
A. Main Results 

The technique of quasi-adiabatic continuation is a general technique for describing 

how the ground state of gapped Hamiltonians evolves under a change in a parameter. We 



consider a parameter dependent Hamiltonian Ti s = J2zhz(s) for some parameter s. We 
assume the ground state or set of ground states is given at s = and we wish to describe 
the evolution of this state or set of states under change in s. We denote the ground states of 
7i s by ^q(s) where a = l...k is an index labeling the k distinct ground states. The ground 
states ^q(s) are not required to be exactly degenerate: they instead each have energy E a (s) 
and we denote the splitting in energy between the highest and lowest energy ground state 
by e. If there is a degeneracy of the ground states at some value of s, then there is some 
freedom in choosing the ground state basis \I/oj we resolve this freedom by requiring that if 
E a (s) = E b (s) for some s, then 

(tt*(s),3,#S(s)) = whenever E a {s) = E b {s), (5) 

where (...,...) is used to denote the inner product. Similarly, we require that 

(*S( a ),0,*s( a )> = o (6) 

for all a, s. We label other states of the system by iff 1 , i > k, with energy Ei(s). 

The most important result of [2] shows that the evolution of the ground state under a 
change in s can be well approximated by a local Hermitian operator acting on the ground 
state if there is a gap. Let there be a gap AE between the sector of ground states and the rest 
of the spectrum for all s. Let ^q{s, sz) denote the set of ground states of the Hamiltonian 
Tl(s) + (hz(sz) — hz(s)). Then it was shown that, for any length lz, there exists a Hermitian 
operator D l ^ c {s) and an ant i- Hermitian matrix N b Z J oc such that D l ^ c {s) has support on the 
set of sites j with dist(Z,j) < lz and satisfies HPjj^s)!! < t q \\d Sz hz(sz)\\, and such that 

^VI>S( S ,^)-^ C (s)^S-HAr^ oc ( s )xI/^( s )| < C 1 (e X p[-Z^/r] + |^|e X p[- > u'Z^])t i? ||^^( S ^) 

b 

(7) 

where C\ is some numerical constant of order unity (throughout, we denote such numerical 
constants by symbols Ci, C2, ■■■) and we define the constants t q by (the q in t q stands for 
"quasi-adiabatic" ) : 

£' = 2v'/AE, (8) 



t q = \jl z /v'AE (9) 
In Eq. ([Tj) and throughout this section, all derivatives with respect to sz- are taken at 



From this, we can define the operator D (s) = XIz^z ( s ) an d bounds on |<9 s \I/q(s) — 
iV loc (s)% - E b N l b ° a c (s)^ b (s)\ follow from Eq. ©, where N l b ° c (s) = Y, z N bl ? oc {s). 

We briefly sketch the derivation of Eq. (JTj) following ideas from [2j. We define the operator 



V z {s) = i / dre-^ 2 / 2 [wl Sz (tr) - h.c] 



(10) 



where w ZjSz = (d Sz h z (s z )) and where for any operator A we define, following 4]: 

i(t) = A(t)exp[-(t/t g ) 2 /2], 
A ± (ir) = — I dtA(t) . 



(11) 
(12) 



The time evolution of operators in Eq. (flT!) is the usual Heisenberg evolution with the 
Hamiltonian 7i(s). Define 



N b Z a (s) = (* b (s),d sz %(s,s z )) - (* b (s),iV z (s)%(s)). 
Then, one may show using the gap that for any \I/g we have 
\d sz %(s, s z ) - iV z {s)%{s) - J2 N b z a (s)V b 

b 
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(13) 



(14) 



E 



Ei(s)-EJs) 



((¥, d sz h z (s z )%))¥ - iV z (s)%(s, s z ) - N b z a ( S )¥ (s) 



< C 2 (t q e W [-(t q AE) 2 /2])\\d sz h z (s z )\\). 

We define T> l ^ c {s) as follows: define Ui to be a unitary operator on the set of sites j such 
that dist(Z, j) > I. Let n(Ui) be the Haar measure on such operators. Define 14 \ 



V l ° c (s) = J dfi(U lz )U} z V z (s)U h 



(15) 



Using Lieb-Robinson bounds, it is possible to show that 



\\V l ° c (s) -T> z (s)\\ < C 3 (exp[-(l z /v't q ) 2 /2] + \Z\exp[-^l z ])t g \\d Sz h z (s z )\\. (16) 

Combining these results (fTHfTBl gives Eq. (ED, where N b z J oc (s) = (^ b (s),d Sz %(s, s z )) - 
(^ b (s),iV z<loc (s)%(s)). Define the operator N z ' loc (s) = £ ba ^° c (s)^) (tfg and N loc (s) = 
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B. Berry Phase 



The discussion that follows in this subsection is meant to introduce the idea of computing 
a Berry phase for quasi-adiabatic evolution and show that it is close to the usual definition 
of the Berry phase. This discussion is not used in the rest of the paper and may be skipped 
if desired, but seems to be interesting in itself as well as for other future applications. 

Define Po(s) to be the projector onto the ground state sector: Po(s) = J2a=i \l/g(s))(\l/o(s). 
A direct computation shows that for all a, b 

\(V b {s),V Ziloc (s)%(s))\ < C 4 (t q ex P [-l z /C]+t q \Z\exp[-fi'l z ]+t 2 q e)s\\d Sz h(s z )\\. (17) 

Note that for large l z and small e, the right-hand side of this equation becomes close to zero, 
giving a result similar to Eqs. ([ofo]) ; this will be used to relate the evolution in the ground 
state sector produced by %D(s) to a Berry phase. 

Eq. ([7]) gives the evolution of wavefunctions under an infinitesimal change in s. We now 
present a unitary operator V to describe the change in wavefunctions from s — to s — 1. 
We define 

V(s) =S'exp[f S ds'tV{s')}, (18) 
Jq 

where S' denotes that the exponential is s'-ordered. Then, ^^(O) is close in norm to some 
linear combination of states ^(s), with an error that may be estimated from Eq. ([7j). 
If the Hamiltonian is quasi-adiabatically continued around a closed path in parameter 
space, so that H(s fmal ) = H(0), then \V(s final )% - E b Qba^ b \ < Ci£z(fflcp[-V£l + 
| Z | exp [-ii'l z ])t q s f ina i || d sz h(s z ) || where 

Q = <S'exp[ f ds'iP (s')T>(s')P (s') - P ^2(d s ,%(s'))) {^(s')] (19) 
Jo 

= S'expf- f ds'N loc (s')}. 
Jo 

Eq. (TT7|) gives 

HQ-Q^II < C^EfoexpHz/f] + t q \Z\ exp[-//Z z ] + t 2 q e)s fmal \\d Sz h(s z )\\, (20) 
z 

where Q berr y is the non-Abelian Berry phase for an adiabatic evolution in a related system 
in which the level splitting between the ground states ^o(s) is set equal to zerojl5]. 

While we used the assumption that there is a gap AE to derive Eq. (p3J, it is important 
to note that this equation in fact depends only on the weaker assumption that the state 
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d s hz{s)^Q is a linear combination of ground states and states ^ l with E{ > AE + E a for 
a = l...k. This follows automatically if there is a gap in the spectrum, but it also follows in 
a spin system if, for example, there is a gap AE between the ground state with spin-0 and 
the lowest spin-1 state, and d s hz(s) has spin-1. 



II. GOLDSTONE'S THEOREM IN SPIN-GAPPED SYSTEMS 



A. Statement of Result and Examples 



Recently, _t 
been proven 



re exponential clustering of correlation functions in a system with a gap has 



n 



m [Si |y|- That is, it has been shown that, for a local Hamiltonian with a 
gap AE between a sector of degenerate ground states, labeled \1/q for a = l...k, and the 
rest of the spectrum, the connected correlated function, (AB) — (AP B), of two operators 
A, B with support on sets X, Y is exponentially small in the distance dist(X, Y). Here, 
(O) = AT 1 Ea=i(*o> °^o) and Po = Eti ^o)(^o is tne projector on the ground state sector. 
For a parameter dependent Hamiltonian, we define (0) s = k~ l J2a=i(^o( s ) , O^o( s ))- in the 
event that the states are not exactly degenerate but only approximately degenerate, 
similar bounds were found with corrections that depend on the energy difference between 
the states. 

In this section we prove a stronger statement about the decay of correlation functions in 
gapped systems with a continuous symmetry. We consider Hamiltonians which obey Eq. (TTJ) 
for some /i, Si and we also require that for any site i the number of sites j within distance / 
of site i is bounded by al d for some finite a, d: 

1 ^ ald - ( 21 ) 

j,dist(i,j)</ 

We consider Hamiltonians which have a U(l) symmetry as follows: for each site i we assume 
that there exists a local operator g« with support on that site such that [Q,Ti] = where 
Q — J2i Qi- F° r an y set X, we define R(9, X) = ILex ex P[^i$]- We consider operators (fix, 4>y 
with support on sets X, Y which transform as vectors as follows under this U(l) symmetry: 
R(-9,X)(fi x R(9,X) = exp[i9)<t> x and R{-9, Y)0 Y R(9, Y) = exp[-i0]0 y . We assume that 
||<?i|| < Qmax for some q m ax- We will show that, for a Hamiltonian which obeys all these 
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conditions and has a gap AE between the ground state sector and the rest of the spectrum 



(<j>x<f> Y ) < C 5 ||</>x0y|||X|adist(X,y) d Sl g ma:E V / dist(X,r)A J E/^x (22) 
(exp[-//distpf, Y)/2] + exp[-dist(X, Y)/2£']). 

Picking the optimal value of // in Eq. we show that 



{<px<P Y ) < C , 6 ||0x0y|||X|dist(X,r) d g m a^ lv /dist(X,r)AE/i;'exp[-dist(X,r)/2^, (23) 
with 

£ = 2v/AE + /i. (24) 

We emphasize that we do not require the states in the ground state sector to be degenerate 
with each other, simply the existence of a gap between that sector and the rest of the 
spectrum. Further, the bounds do not depend on the value of k or on energy difference 

nn 

between the states This result ( 1221(231) is stronger than that in [5|, ll6( which was only 
valid for d < 2 while the present result is valid in arbitrary dimension. It is also stronger 



than other previous results 17j which either required a unique ground state or else assumed 
an ergodic property which is equivalent to requiring the vanishing of the matrix elements in 
the ground state sector; in either of these cases the decay of correlations becomes equivalent 
to clustering. 

Before giving the proof we motivate the definition of qi by discussing physical examples. 
In a Bose system with conserved particle number, the can represent the particle number 
on a given site and the operators 0x,0y can represent creation and annihilation operators 
for the bosons. For a spin system, the can represent the z component of the spin on a site 
and the 0x,0y can represent raising and lowering spin operators on sites. Consider such 
a spin system in a disordered quantum paramagnetic phase. An example of such a system 
would be obtained by a two dimensional Hamiltonian of spin-1/2 spins for a system of two 
layers with coordinates labeled (i, I) where / = 1, 2 is a layer label and i indexes position in 
each layer. First consider the Hamiltonian: 

H = J S it i' S jt i + J± ^2 S itl ■ S ij2 , (25) 

i,j n.n. i 

with J, Jj_ > where the first sum is over nearest neighbor This is a bilayer model 



Heisenberg model that is connected with the Kondo lattice model 18|. For J± » J, the 



S 



FIG. 1: Illustration of the geometry we consider. X, Y are shown as shaded regions, while X' 
includes everything within the outer circle around X. 



ground state is unique and gapped, with spins (i, l),(z, 2) being in a singlet state with 
high probability. In this case, Eq. ( 122]) shows (Sf^Sj^) is exponentially decaying in dist(z, j) 
but this provides no additional information beyond that already known from the exponential 
decay of correlation functions. Now consider the Hamiltonian on a lattice with two "defects" : 
for some given k, I, we replace the spins at sites (k, 2) and (Z, 2) with spin-0 spins 19J). That 
is, we remove them from the lattice, so Sk,2 = Sip = 0. Then, there are unpaired spins at 
(k, 1) and (1, 1). In the limit J = 0, indeed there are four exactly degenerate ground states. 
Returning to the case J± » J > 0, we expect that if dist(/c, I) is large then there will be 
four low energy states and then a gap to the rest of the spectrum. In this case, Eq. ( 122]) 
can be used, given the assumption of a gap as all other requirements are trivially satisfied, 
to bound (SkiSfi) and hence while the operators S^i,^,! may have nonvanishing matrix 
elements between these states, the average of the correlation function over the different 
ground states is small 



B. Proof 

To show Eq. ( 1221 . we define a set of parameter dependent Hamiltonians He = J2z hz(9) 
as follows. Let X' denote the set of sites i such that dist(X, i) < dist(JT, Y)/2, as shown in 
the figure. Then define h z (0) = R(X',6)h z R(X',-6). Clearly, as R(X',-9) is a unitary 
transformation, He has the same spectrum of H and the ground states of He are given by 

%(6) = R(X',9)M0). 
Thus, 

de(4> x Y )e = de{R{X',-e)4, x ^ Y R{X',e)) (26) 
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= d e exp[i9]((j)x(pY) = *(0jc<M» 

where we used the fact that X G X' while Y D X' = so that [0 y , R(X', 9)] = and where 
we evaluate the derivatives at 9 = 0. 

Without loss of generality we may assume that 

[Q,h z ]=Q (27) 

for all Z. To show this, we use [Q,T~i\ = to write 7i = J2zOz where Oz = 
(l/27r) Jq 77 d9 exp[iQ9]hz exp[— iQ9]. Then, Oz has support on Z and ||O z || < \\hz\\- Fi- 
nally, [Q, Oz] = 0. Hence, by replacing hz by the operators Oz, we succeed in rewriting 
the Hamiltonian such that Eq. (127]) is satisfied and such that the Hamiltonian still obeys 
Eq. p. From Eq. (1271) . 

Znl' = -> fcs(0) = /i z , (28) 
Z C X' -> fc z (0) = /i z . 

That is, hz(9) = hz unless Z contains some point i E X' and some other point j ^ X'. 

We now use quasi-adiabatic continuation to bound the first line of Eq. (126|) . thus 
bounding the correlation function on the last line of this equation. Using dg(<j)x4>Y)d — 
lEa = ide(%(8)Ax0Yn(Q)} and Eq. © we have 

k 



On 



)o\ < |ETEOWWJ0x0y,^ = O)]^))+ (29) 



Z k a: 1 



£<' Zoc (0 = o)(*g(e),^*S(e)> + Y,^ loc {9 = o)*(*S(e),^^*S(fl)> 

6=1 6=1 

+C 1 ]T(exp[-/ z /£ / ] + |Z|exp[-/i7z])^ll^z(0z)|| 
z 

The terms involving the matrix N z,loc in Eq. (1291) cancel after summing over a and b. 
This is the key step which makes the quasi-adiabatic continuation useful in deriving the 
desired exponential decay of correlations. We will now show that all the other terms on 
the right-hand side of Eq. (1291) are exponentially small in distpf, Y), when we choose 
lz = min(dist(X, Z), dist(Y, Z)) — 1. This requires a straightforward, but lengthy, series 
of inequalities. 

The commutator T) l ^ c {9 = 0)] is vanishing unless T>z c (9 = 0) has support on either 

set X or set Y which, with the above choice of lz, occurs only ifZflX^0orZny^0. 
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Further, from Eq. (|28|) V l ^ c (6 = 0) is vanishing unless Z includes some point i G X' and some 
other point j G" X'. Thus, the commutator [</>x</y, T> l ^ c (6 = 0)] is vanishing unless either Z 
includes some point i G X and some other point j G" X' or else Z includes some point i E X' 
and some other point j G Y. In either case, we have i G X' and dist(i, j) > dist(X, Y)/2. 
Thus, 

EII[^.^ = o)]ll < ^UxM E E E (so) 

2 ieX' j,dist(ij)>dist(X,y)/2 

< 4t,g maa; ||0 x y |||X / | Sl exp[V(dist(X,F)/2)] 

< 4t g g max ||0 x y ||a(dist(X,r)/2) d |X| Sl exp[- yU '(dist(X,r)/2)], 

where we used Eqs. JTEH) and ||D| C || < t ? ||0* z /i z (0z) || and ||d^/i z (0 z )|| < 2q max \\h z \\ for 
the last two inequalities. 

Similarly, noting that l z > dist(X, Y)/2 — diam(Z), 



z 

< 

iex' Z3i 



E|(exp[-/z/r] + \Z\exp[-fj,'l z ])\t q \\d ez h z (O z )\\ (31) 
z 

E E 2(exp[-(dist(X,r)/2 - diam(Z))/£'] + |Z| exp[-//(dist(X, F)/2 - diam(Z)])t, gmaa! ||/i z || 



< \X\a(dist(X,Y)/2) d s 1 q max ^dist(X,Y)AE/v'(ex.p[-fi'dist(X,Y)/2] + exp[-dist(X, y)/2£']J. 
Combining Eqs. fl30|31l) . we find that 



\do(<f>x<h)e\ < C 5 U x 4> Y \\\X\adist(X,Y) d Sl q max ^dist(X,Y)AE/v> x (32) 
(exp[-//dist(X, Y)/2] + exp[-dist(X, K)/2£']). 

Combining this with Eq. (|26|) . we arrive at Eq. (l22l) . as desired. 

It is interesting to note that Eq. (|2"T1) is necessary in this derivation. We sketch a system 
for which Eq. (12 lj) does not hold, and show how a Goldstone theorem may fail in this case. 
Consider a random graph with V nodes each having coordination number 3. Consider a set 
of V spin- 1/2 spins, with Hamiltonian 

H — — E J'y'S'i ■ Sj, (33) 

where the interaction matrix equals 1 if the nodes i,j are connected by an edge on the 
graph, and zero otherwise. The interaction is ferromagnetic, so pointing all spins up (or 
in any other direction) gives a ground state. Further, the Hamiltonian is local, using a 
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shortest path metric on the graph to define dist (z , j). However, a random graph of this form 



is typically an expander graph 20| with a gap in the spectrum of the graph Laplacian, so a 
spin-wave theory calculation 2l| gives a gap in the magnon spectrum. Thus, this system has 
a set of degenerate ground states and a gap. However, the spin correlations do not decay, 
as (Si ■ Sj) = 1/4 for all i ^ j. 

III. FLUX QUANTIZATION IN SPIN-GAPPED SYSTEMS 
A. Examples and Statement of Result 

In this section we consider flux periodicity in systems of fermions with half-integer spin. 
Before stating the precise result we motivate some of the definitions through physical exam- 
ples. Consider a system of fermions on a cylinder, and suppose that an additional magnetic 
field is applied through the cylinder, introducing an Aharonov-Bohm phase for the fermions. 
The Hamiltonian is unchanged, up to a gauge transformation, if the magnetic flux through 
the cylinder is changed by 2ir, and hence the ground state energy is periodic in the flux with 
periodic 2tt. For an arbitrary Hamiltonian the energy of the ground state for tt flux may be 
much different for that at flux. However, a superconductor typically has its ground state 
energy approximately periodic in the flux with period it due to so-called Byer-Yang states s| , 
up to some corrections which are exponentially small in system size L. This is often taken 
to imply that the elementary excitation in the system has charge 2e rather than e. 

However, there are other interesting scenarios in which the ground state energy at 7r flux 
is close to that at zero flux. Consider a system at half filling which forms a resonating valence 
bond state (RVB) with a spin gap. Then dope this system by introducing a small number of 
holes. The holes have vanishing spin and charge +e, and may Bose condense, again giving 
a superconductor. Naively, one may expect that in this case the ground state energy of 
the system at 7r flux would be very different from that at zero flux as the excitations that 
condense have charge e. However, the presence of topological excitations of the valence bond 
system implies that in fact the ground state energy is still approximately periodic in the flux 
with period 7rj^, Q, 23]. The question was raised in whether there were some general 



conditions under which one could prove that this approximate periodicity held for fermionic 
system. In this section, we provide a partial answer to this question by considering the case 
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in which the system has a spin gap, as would occur for an s-wave superconductor or for a 
short wave RVB state, and prove that there exists a state of the system at tt flux which is 
close to the ground state energy at zero flux. Unfortunately, our results will not apply to 
d-w&ve superconductors where there are gapless spin excitations. 

We now state the result. We assume that there are conserved spin and charge as follows. 
On each site i we assume that there is a Hilbert with dimension 4, representing the states 
of no fermions, one fermion with spin up, one fermion with spin down, and two fermions. 
We define the spin operator on site i, S?, to be diagonal in this basis with eigenvalues 
0, +1/2, —1/2, respectively, and we define the fermion number operator on site i, qi, to 
be diagonal in this basis with eigenvalues 0,1,1,2. We assume that the Hamiltonian 7i 
commutes with the operators S z = J2i Sf and Q = J2iQu as before, this means that we may 
assume without loss of generality that [S z , hz] = [Q, hz] = for all Z. Note that Sf ± qi/2 
has integer eigenvalues. We assume that the Hamiltonian is local in the sense of Eq. (pQ) 
and we assume that the system is periodic in one direction as follows: let the system be 
defined on a d dimensional lattice of V different sites, with any site i having coordinates 
(xo(i),xi(i), ...Xd-i(i)), with periodic boundary conditions in the xo direction with period L. 
Assume that dist(i, j) > min n \xo(i) — Xo(j) — nL\, where the minimum is over all integers 
n. This is satisfied by any physically reasonable metric on the lattice, such as a Manhattan 
metric. 

We assume that TC has a unique ground state although the results can be readily 
extended to a system with multiple ground states with small splitting e. We define a Hamil- 
tonian 7i to have a spin gap AE if the the ground state has vanishing spin and the first 
excited state with non-vanishing spin has energy at least AE above the ground state energy. 

We now define the Hamiltonians Tie with flux 6 inserted. Define £ to be the set of 
all points i with < xo(i) < L/2. For any set X, define R q (9,X) = Tiiex exp[iqi9]. 
Define Tig = Y^zhz{@), where we define hz as follows: if Z contains any points i with 
L/A < Xo(i) < 3L/4, then we define 

h z (9) = R g (-9)h z R g (9) (34) 

Otherwise, we define 

h z {9) = h z . (35) 

We will prove that, under the conditions identified above, (conserved spin and charge, 
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x=0 x=L 

FIG. 2: Illustration of the geometry we consider. Flux is inserted along the dashed line in the 
center, and at the edge. The slanting upwards line is the support of Wi, the slanting downwards is 
the support of Wi- The grey areas have width L/2(/z£+l). The lines xq(z) = L/A and xq(i) = 3L/4 
pass vertically through the middle of the left and right grey areas, respectively. 

half-integer spin, locality, periodicity, uniqueness of ground state, and spin gap), that there 
exists some state, which we will write as W^ir^o with W\ a unitary operator, such that 
the difference in energies 

KWi(7r)tt ,W*Wi(7r)tt ) - (*o, < C 7 V 2 y/L/vAEs* exp[-L/4(£ + 1///)]. (36) 

Note that Wi(ti)^o need not be the ground state of Ti^. 



B. Proof 

In this subsection we provide the proof. The main idea is to define separate fluxes for 
up and down fermions separate, and use quasi-adiabatic continuation to carry one flux from 
to it and the other from to —71. Then, the proof closely follows the proof of a higher 



dimensional Lieb-Schultz-Mattis system given in 24]. For any set X, define R^(6,X) = 
Eliex eyL~p[i{Sf + qi/2)9] andi?j(^,X) = Y[ i&x exp[i(— S!f +g»/2)0]. Thus, produces a gauge 
transformation on the up spin fermions and Ri transforms the down spin fermions. Note 
that R^(tt, X)Ri(— it, X) = R g (ir,X); this equality depends on the fact that 7r, X) = 
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We now define a family of Hamiltonian with separate up and down fluxes inserted, and 
we also insert flux at two different points: along a line with xq = and along a line with 
Xq = L/2 as shown in the figure. By gauge invariance the spectrum of the Hamiltonian 
depends only on the total of these fluxes, but the introduction of two different flux angles is 
a useful technical device. We define four flux angles, &\, &[, &l B 2 , where the fluxes described 
by angles B\ are along the line with xq = and the fluxes described by angles 02 are along 
the line with xq = L/2, and where we introduce different fluxes for up and down particles. 

I T 

The insertion of physical magnetic flux into into the system the system leads to 9{ = 9[ and 
#2 = #2, however our introduction of separate angles for up and down particles is necessary 
to derive the final result. The resulting Hamiltonian H(B\, B[, B\, B l 2 ) = £ z h z (0\, B[, Bl B l 2 ), 
where we define h z as follows: if Z contains any points i with L/4 < Xq(i) < 3L/4, then we 
define 

h z {9\, B{, el 9i) = R^(-9l C)R L (-0i, C)h z Ri(9l, C). (37) 

Otherwise, we define 

hz{el 9{, el el) = R^el QRtfi, c)h z Rr(-el c^i-el c). (38) 

Then, if e\ = -el and d{ = -dl we have ^0(^1,^,-^,-^1) = R^ifil, C)Ri(Pi, ^)^o- 
Further, the eigenvalues of the Hamiltonian TC(9\, el Bl el) depend only on e\ + Q\ and 
e\ + el We define Z c to be the set of Z such that h z {B\, e{, e\, B 2 ) is different from h z . 
Crucially, the physical Hamiltonian with 7r magnetic flux inserted is 

H n = H(n,n,0,0). (39) 



The proof now closely follows 24j, and we only sketch the proof; the only differences 



from [24| are the presence of separate flux angles for up and down fermions, and some minor 
complications due to exponentially decaying interactions, rather than finite range. We define 
Wi(4>) to be the operator which quasi-adiabatically continues T~C(9, —0,0,0) from B = to 
= 0, using Eq. (JE) with l z = L/4 - diam(Z) - L/4(£fi + 1). We define W 2 {<j>) to be the 
operator which quasi-adiabatically continues 7i(0, 0, B, —B) from B = to B = —(p. We define 
W{4>) to be the operator which quasi-adiabatically continues H(B, —B, —B, B) from B = to 
= 0. 

n 

Following the results in [24J one may show that 

\R^,X)R l (-7r,X)V - W(rr)^ \ < c 2 , (40) 
15 



where 



c 2 = C s £ t q \\h z \\(exp[-l z /AE/2v'] + \Z\exp[-n'l z )) (41) 

z,zez c 

< Cs^L/vAE Sl V(exp[-l z /AE/2v'} + \Z\ exp[-fi'l z ]). 

and also that 

\\Wi(ir)W 2 (-ir)-W(7r)\\<c 3 , (42) 

where 

c 3 = C 9 J2 t q \\h z \\(&v[-l z /AE/2v^ + \Z\exp[-ti'l z ]). (43) 

z,zez c 

Here, we use constants c 2 , C3 to follow the notation of (2^]. 

Eq. (I42p relies on the fact that Wi(ir) and can be written as exponentials of local 

operators T>i,T> 2 . The commutator of T>\ with D 2 can be bounded, and thus one can "re- 
order" the quasi-adiabatic evolution so that instead of first evolving with W 2 and then with 
Wi, both flux angles are changed at the same time, as with W. 

Eq. (l4"0j) relies on the existence of a spin gap in the system. For a gapped sys- 
tem with a unique ground state, the quasi-adiabatic evolution approximately evolves the 
ground state of the initial Hamiltonian, into the ground state of the final Hamilto- 
nian, R^(tt, X)R±(—tt, X)^ . Here, we have not specified that the Hamiltonian has a gap. 
However, dgTC(9, —9, —9, 9) is a spin-1 operator, and hence the spin gap suffices to show 
Eq. (HO]). 

We now show that (W 1 (n)^ ,'H( y Tc,7i,0,0)W 1 (7r)^ ) is close to (4> ,7i4> ), which is the 
desired result as by Eq. (1391) it shows that W\{tt)^q is a state for which the expectation 
value of the energy with Hamiltonian 7i n is close to the expectation value of the energy of 
the state \l/o with Hamiltonian 7i(0, 0, 0, 0). 

To show that the expectation value of the energy is close, we write 
?i(vr, Ti, 0, 0) = J2z h z (n, ir, 0, 0). We define Si to be the set of points % with 
-L/4 + L/4(^/x + l) < x (i) < L/4-L/4(£/i + l), and define S 2 to be the set of points % with 
L/4 + L/4(£/i + l) < x (i) < 3L/4-L/4(f// + l). Note that W 1 is supported on St and W 2 is 
supported on S 2 as shown in Fig. 2. li ZnS 1 = 0, then (Wi(7r)* , h z (ir, it, 0, 0)Wi(7r)* ) = 
(*o,Wo>. If ^ n 5 2 = 0, then |(Wi(7r)* , M* - . tt, 0, 0)Wi(7r)tt ) 



16 



KWi(tt) W 2 (7r)* , h z (n, vr, 0, 0) Wi(tt) W 2 (7r)tf > 



(^( 7 r,X)^(-7r ) X)* ,^(7r,7r,0 ) 0)it:T( 7 r,X) J RH-7r,X)*o}| < 2(c 2 + c 3 )\\h z \\. Fi- 



nally, the sum of \\h z \\ over all Z such that Z fl 5i 7^ and Z fl 5 2 7^ is bounded by 
Vsi exp[— /iL/4(^yU + 1)]. Thus, we find that 



giving Eq. (136|) as claimed. 

This calculation can be extended to systems with larger symmetry groups than the U(l) 
symmetry above. If a system has N species of fermions, with an SU (N) symmetry, one may 
define angles 01,9% for a = 1...N for each species of fermion. We perform a continuation 
with 91 = 9/(N-l) for 1 < a < N-l, and 0f = -0, going from = Oto0 = 2n(N-l)/N. 
Then, after the continuation we have 01 = 2tt/N for 1 < a < N, and 0± = 2ir/N — 2tt. 
Thus, if the ground state of the Hamiltonian with no flux is an SU(N) singlet and the system 
has a gap to the lowest state which is not an SU(N) singlet, then there is a state of the 
Hamiltonian with flux 2n/N which is close in energy to the ground state of the Hamiltonian 
without flux. Colloquially, the flux is quantized in units of 2n/N. 

IV. DISCUSSION 

We have used quasi-adiabatic continuation as a tool to study systems with continuous 
symmetries and a gap. This leads to a version of Goldstone's theorem for nonrelativistic 
systems which is valid in arbitrary finite dimension and which does not depend on additional 
assumptions regarding ergodicity. While this result was derived from an infinitesimal quasi- 
adiabatic continuation, the continuation under a non-infinitesimal change in parameters 
was used to derive an approximate flux periodicity of systems with a spin gap. While this 
periodicity is realized in different ways in different systems, for example by Byers-Yang 
states in an s-wave superconductor and by topological excitations in an RVB state, the 
mathematical result covers both cases. 

Finally, we have discussed the Berry phase under quasi-adiabatic continuation, and shown 
that it is close to the usual non-Abelian Berry phase for adiabatic evolution. This result 



I (Wi(7r)*o, H(tt, it, 0, 0)Wi(7r)*o) - (*o,W* >l 



(44) 



< C 7 V 2 ^L/vAEslexp[-(L - L/Afo + !))/£] 



= C 7 V 2 yjL/vAEsl(exp[-L/4(S + V^)}), 
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will be used in a future work to consider Hall conductance quantization in many-body 
systems, and, under the assumption of an excitation gap, to remove the need for an averaging 
assumption on the Hall conductance 25]. 
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